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ABSTRACT 

We present results of spectral analysis of Ginga data obtained during the decline 
phase after the 1989 outburst of GS 2023+338 (V404 Cyg). Our analysis includes 
detailed modelling of the effects of X-ray reflection/reprocessing. We have found that (1) 
the contribution of the reprocessed component (both continuum and line) corresponds 
to the solid angle of the reprocessor as seen from the X-ray source of Q ss (0.4 — 0.5) x 2ir, 
(2) the reprocessed component (both line and continuum) is broadened ("smeared") by 
kinematic and relativistic effects, as expected from the accretion disk reflection. We 
discuss the constraints these results give on various possible system geometries. 

Subject headings: accretion, accretion disk - black hole physics - binaries: close - stars: 
individual (V404 Cyg) - X-ray: stars 

1. Introduction 

Some of the strongest evidence for the existence of accretion disks around black holes has 
come from X-ray observations of the relativistically smeared iron Ka line profile in Active Galactic 
Nuclei (Tanaka et al. 1995; Iwasawa et al. 1996; Nandra et al. 1997). This fluorescence line is 
produced by hard X-ray illumination of the accreting material, and the combination of high orbital 
velocities and strong gravity in the vicinity of a black hole gives the line a characteristically skewed, 
broad profile (Fabian et al. 1989; Laor et al. 1991). A reflected continuum should also accompany 
this line (e.g. Lightman & White 1988; George & Fabian 1991; Matt, Perola & Piro 1991), 
and the amplitude of both reprocessed components gives constraints on the solid angle subtended 
by the accreting material, its inclination, elemental abundance and ionization state. Both the 
amount of reflection/fluorescence and the intensity of the relativistic effects on the observed line 
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profile strongly support the idea that the accretion disk extends down to the last stable orbit in 
AGN. The Black Hole Candidates (BHC) show many X-ray spectral similarities to AGN, plausibly 
because both involve the same physical processes of disk accretion onto a black hole (Inoue 1993; 
Ueda, Ebisawa & Done 1994; Gierlihski et al. 1997). 

Many Soft X-ray Transients (SXT) are known to be BHC. Their mass accretion rate varies 
over at least 2 orders of magnitude from outburst (where it is at about Eddington) through the 
decline on time-scales of months (see Tanaka k, Shibazaki 1996 for recent review). A number of 
these systems were observed by Ginga and the spectra obtained are amongst the best in which to 
investigate the overall effects of X-ray reprocessing. Moderate spectral resolution (18% at 6 keV) 
is more than compensated by very high signal-to-noise, broad band pass (from 1 keV up to 20-30 
keV) and our relatively good understanding of the instrument. 

In this Letter we present results of spectral analysis of data obtained during the 1989 outburst 
of GS 2023+338 (V404 Cyg; Kitamoto et al. 1989). The outburst was well covered by Ginga from 
its initial detection by All Sky Monitor on May 22 until November 1989 (Tanaka & Lewin 1995). We 
have selected two data sets where there is little short time-scale spectral variability, obtained June 20 
and July 19-20, respectively. We analyse these using models of X-ray reprocessing that consistently 
connect the properties of the iron Ka line with the properties of the reflected continuum. We 
show that the reflected component is present in the spectra, that it is smeared by kinematic and 
relativistic effects as expected from an accretion disk, and that both its normalization and amount 
of smearing imply that the disk does not extend down to the last stable orbit. 

A full analysis of the data covering the outburst and entire decline phase will be presented in 
future paper (Zycki et al., in preparation, hereafter paper II). 

2. Data selection and reduction 

The data were extracted from the original FRF's using the Ginga reduction software at Le- 
icester University. Paradoxically, background subtraction poses a problem for a source as bright as 
GS 2023+338 as the background monitors are strongly contaminated by source counts. However, 
the background can be estimated from nearby observations at similar points in the satellite orbit, 
and its fractional contribution is low (< 3% below 10 keV). Full details of our method of background 
subtraction will be given in paper II. We allow for 0.5% systematic error in the data. 

3. Model 

The model components are: multi-temperature accretion disk spectrum ("disk black body"; 
Mitsuda et al. 1984) and a power law with its Compton reflection from a possibly ionized medium 
including the iron Ka line emission. The Compton reflected component is computed using the 
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XSPEC version 9.01 model "pexriv" (Magdziarz & Zdziarski 1995) with the ionization parame- 
terised by the ionization parameter, £ = Lx/nr 2 as in Done et al. (1992). We assume elemental 
abundances of Morrison & McCammon (1983) with the exception of the iron abundance, [Fe], 
which is a free parameter. The iron Ka line is computed using modified Monte Carlo simulations 
code of Zycki & Czerny (1994). We updated the Fe Ka line energies, fluorescent yields and Fe 
K-shell edge energies after Kaastra & Mewe (1993). 

The reprocessed component can then be "smeared" to simulate the relativistic and kinematic 
effects of disk emission (see e.g. Fabian et al. 1989; Ross, Fabian & Brandt 1996). We assume 
a non-rotating black hole so the model is parameterised by the inner and outer radius of the disk, 
R m and .Rout respectively, and the radial distribution of irradiation emissivity. We fix the form of 
the emissivity, F 1TI (r) cx r~ 3 as expected from coronal illumination, fix i? Q ut a t 10 4 R g , since it is an 
appropriate value for LMXRB (R g = GM/c 2 ) and fit i?j n and the normalization of the reprocessed 
component, / = £l/2ir, where fi is the solid angle of the reprocessor as seen from the X-ray source. 
We emphasize that i?i n , R out and Fi rr (r) describe only the smearing effect, not CI. 

We assume that the inclination of the system is i = 56° (Pavlenko et al. 1996). 

4. Results of model fitting 

We begin with the simplest possibility, that is a power law spectrum with a narrow line at 6.4 
keV. We than add the reflected continuum assuming first that it is un-ionized, non-smeared and 
the iron abundance is 1 relative to cosmic value. In the next step we fit £ and [Fe] and finally we 
introduce the effect of smearing. This multi-step procedure is important (1) to demonstrate the 
significance of reprocessing and (2) since the moderate spectral resolution of Ginga means that the 
broadening of iron spectral features due to ionization can mimic the relativistic smearing. 

For the June 20th data set, the simplest model (model in Table 1) gives = 1.9 and it is 
thus not acceptable. Adding the reprocessed component is highly significant even in the simplest 
version (model A), \ 2 = 19.9/25 d.o.f. (x 2 = 0.80). We note that the normalization of the reflected 
component is significantly smaller than 1, / = 0.46 =L 0.04. 

The fit can be improved by allowing for ionized reflection (model B). Assuming the reprocessor 
temperature T = 10 6 K (for the purpose of computing ionization balance only), the best fit has 
X 2 = 15.2/24 d.o.f. for £ = 0.1 +4 . If, instead, we let iron abundance be free whilst fixing £ = 
(model C), we obtain x 2 = 19.1/24d.o.f. for [Fe] = 1.20±0.35, i.e. [Fe] is not a significant parameter. 
The best fit with both £ and [Fe] free (model D) has x 2 = 14.9/23 d.o.f. for £ = 0.1 but again [Fe] 
is consistent with the cosmic value. 

We now introduce the effect of smearing (model E; both £ and [Fe] are left free as well). This 
results in a decrease of x 2 by Ax 2 = 8.9 which is significant at more than 99.9% confidence level 
(the F-test for one additional parameter). The best fit inner radius is R m = 25±ff i? g . Figure | 
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shows the spectrum and residuals of the best fit models D and E. 

A similar progression in quality of the fit is given by the July spectrum (Table 2; Figure 1), 
where again the relativistic smearing effects are significantly present in the data. 

We have also tried a model with £ changing as a function of radius, £(r) oc r a , but with no 
relativistic smearing. Best fit of the model for the June data has \ 2 = 14.6/22 d.o.f., while for the 
July data \ 2 = 20.1/22 d.o.f . 

5. Discussion 

5.1. Primary continuum 

From fits we have T ~ 1.7, while contemporaneous high energy data show a cutoff at about 100 
keV (Sunyaev et al. 1991). The source must then be marginally optically thin if the primary X-rays 
are produced by thermal Comptonization, with tt ~ 0.5 and 1.5, in a disk or sphere geometry, 
respectively (Titarchuk 1994). 

5.2. Geometry 

The results of spectral modelling clearly show that the reprocessed component is present in 
these data. The required smearing cannot be explained by a radial distribution of ionization 
since the data strongly favour a low and uniform ionization. Thus the relativistic effects in an 
accretion disk/black hole system are the most plausible explanation. The reflected fraction for 
both continuum and line is roughly ~ 0.5 x that expected from an isotropically illuminated flat 
disk. Thus the line is not depleted by Auger ionization masked by relativistic smearing (Ross et 
al. 1996). The line is weak because the covering fraction of the reflecting material is small, and 
Auger ionization is ruled out by the low ionization state of the disk. The low covering fraction does 
not seem to be an artifact of super-solar abundances: fixing / = 1 and allowing the abundance 
to be free results in a poorer fit, with x 2 = 19/23 d.o.f. and 16/23 d.o.f. for the June and July 
data, respectively, and the overall abundances, A ~ 8. It is also commonly seen in the persistent 
BHC: Cyg X-l and GX 339-4 (Done et al. 1992; Gierlihski et al. 1997; Ueda et al. 1994), so 
is not some time dependent effect of disk evolution in transient systems, but rather represents 
a significant geometrical constraint. Another constraint comes from the fact that the observed 
relativistic smearing is less than that expected from a central point source illuminating a flat disk 
which extends down to 6R g . 

One possible geometry is a spherical or flattened X-ray source centered on the black hole 
as originally proposed by Thorne & Price (1975) for Cyg X-l. A spherical source with tt 3> 1 
gives / = 0.5 naturally (Done et al. 1992), as the X-ray emission only escapes into a hemisphere 
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tangential to the source surface. However, here the source is probably marginally optically thin 
with tt ~ 0.5 — 2 (see § ^^). This can still give a reduction in / since there is a hole in the inner disk 
at r < Ri n . Some of the photons escape without illuminating the disk, and the remaining reflection 
fraction is diluted by hard X-ray photons produced on the other side of the disk. The relativistic 
smearing constraints then imply that the the source emissivity is less steep than oc r~ 3 and/or the 
inner disk disk truncates at R{ n > 6R g and/or the outer disk flares. There are few real physical 
constraints on the source emissivity. Local release of gravitational energy in a Keplerian disk 
gives a luminosity oc r -3 , but non-local mechanisms could also operate, perhaps even producing 
a constant source oc r°. Outside of the source the illumination is approximated by a central point 
source i.e. oc r -3 but within the spherical source the illumination produces an emissivity of oc r~ 2 
(for luminosity oc r~ 3 ). This can reduce the most strongly smeared components to the observed 
level, though truncation of the inner disk radius may also be required. The geometry might then 
be in accord with recently proposed scenario of SXT evolution after outbursts (Esin, McClintock 
&: Narayan 1997 and references therein) although the overall behavior of the source does not seem 
to follow the proposed scenario. 

Small active regions on the disk (a "patchy" corona), perhaps powered by magnetic recon- 
nection give another possible geometry (Haardt, Maraschi & Ghisellini 1994, Stern et al. 1995). 
Physical support for such a picture may come from recent models of accretion disk viscosity as 
an MHD disk dynamo (Hawley, Gammie & Balbus 1996). However, these models give / ~ 1 
unless the height of the reconnecting regions, h rec , is large compared to the size of the disk. This 
would require that the reconnecting regions be concentrated along the inner disk radius R[ n , and 
that h rec ~ i?j n , since the outer disk is expected to be very large. This would also go some way 
to satisfying the smearing constraints, as the irradiation would then be constant from ~ i?; n to 
~ i?i n + h rec rather than oc r~ 3 . 

Another geometry that has been proposed, though with rather less physical motivation, is 
a continuous corona overlying the accretion disk (Haardt & Maraschi 1993). A fraction ~ 1 — 
exp(— t I cosi) ~ 0.5 of the reflected spectrum is itself Compton scattered by the X-ray emitting 
corona, and loses its characteristic spectral shape (Haardt et al. 1993). The energy generation 
would be expected to be oc r~ 3 , so the smearing constraints are a problem unless the inner disk 
truncates, or the outer disk flares. However, this simple scenario can be ruled out from detailed 
spectral fitting of Cyg X-l (Gierlihski et al. 1997; Poutanen, Krolik & Ryde 1997), although perhaps 
this merely points to added complexity such as temperature/optical depth structure in the corona. 

Additional independent constraints on the geometry are provided by low normalization of the 
soft component and weak ionization of the reflecting medium. They both support the idea of the 
disk being truncated at ~ 30 R g as then both the expected temperature of Shakura-Sunyaev (1973) 
disk emission component (for m = 0.01 MEdd an d M = 10 Mq), T ~ 0.15 keV and the ionization 
parameter, £ ~ 1, would be in agreement with our results. However, low £ may also mean that the 
disk is much denser than the SS solution due to e.g. coronal dissipation of energy (e.g. Svensson & 
Zdziarski 1994). 
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6. Conclusions 

We have performed spectral analysis of Ginga data of soft X-ray transient GS 2023+338 
(V404 Cyg) obtained one and two months after its 1989 outburst. We have found that 

• the Compton reflected continuum and iron fluorescent Ka line are present in the spectrum, 

• the properties of both reprocessed components (normalizations, ionization parameters) are in 
agreement, 

• the data require both reprocessed components to be broadened and smeared, 

• the smearing is consistent with being due to reflection from a Keplerian disk, 

• both the normalization of the reflected continuum and the amount of smearing constrain the 
geometry. 

This research made use of data obtained from the Leicester Database and Archive Service at 
the Department of Physics and Astronomy, Leicester University. CD acknowledges support from 
a PPARC Advanced Fellowship. Work of PTZ was partly supported by grant No. 2P03D00410 of 
the Polish State Committee for Scientific Research. 
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Fig. 1. — Best fit spectra (models E in Tables 1. and 2.; upper panels), their x 2 residuals (middle 
panels) and the \ 2 residuals for models D i.e. without kinematic smearing. The dotted histograms 
show the background. Left panels (a) are for June data, right panels (b) for July data. 
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Table 1. June 20th data 



model 


fcTkeV 


r 




[Fe] 


/ Rin (Rg) 


EW (eV) 


X 2 /d-o.f. 



A 
B 
C 
D 
E 


1.51 ±0.08 
0.36 ±0.04 
0.36tg;» 

o.38lS:S 

37+° 18 

u -°' -0.05 

n SQ+ 0,25 

u ""-0.CS 


1.34 ±0.02 
1.66 ±0.01 

1.64 ±0.04 

1.65 ±0.03 
1.64 ±0.04 


0(f) 
0.1+ 4 

0(f) 
0.1+ 4 

+5 


1(f) 
1(f) 
1.2 ±0.35 
1 0+ ' 4 

1 - u -0.2 

1 5+ - 7 

i -°-0.8 


0.46 ± 0.04 

n 45+0 4 

0.43 ± 0.06 
0.44 ± 0.06 
0.47 ± 0.06 25ti\ 


90t?2 


46.8/25 
19.9/25 
15.2/24 
19.1/24 
14.9/23 
6.0/22 
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Table 2. July data 



model kTkeV T £ [Fe] a / R in {R e ) EW (eV) x 2 /d.o.f. 



1.4 ±0.1 1.45 ±0.25 ••• ••• ••• ••• 93 ± 14 57/25 

A 0.36 ±0.03 1.72 ±0.015 0(f) 1(f) 0.40 ± 0.05 ••• ••• 34.6/25 

B 0.441°;*° 1.68 ±0.04 13+f 2 1(f) 0.32 ± 0.08 ••• ••• 22.9/24 

C 0A2+°^° 1.68 ±0.03 0(f) 1.71qJ 0.35 ± 0.05 ••• ••• 26.9/24 

D 0.42±oo6 1 - G8 -oo6 °- 1+4 ° 1-5 ±0.6 0.35±g?g ••• ••• 21.5/23 

E 0.47tg : gjj 1.67t° : °* 0+ 20 2±£ 8 0.37t° : ^ 35±l{° ••• 15.4/22 



a Upper limit on [Fe] is 3 



